Dynamics of single neutral impurity atoms immersed in an ultracold gas 
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We report on controlled doping of an ultracold Rb gas with single neutral Cs impurity atoms. 
Elastic two-body collisions lead to a rapid thermalization of the impurity inside the Rb gas, repre- 
senting the first realization of an ultracold gas doped with a precisely known number of impurity 
atoms interacting via s-wave collisions. Inelastic interactions are restricted to a single three-body 
recombination channel in a highly controlled and pure setting, which allows to determine the Rb- 
Rb-Cs three-body loss rate with unprecedented precision. Our results pave the way for a coherently 
interacting hybrid system of individually controllable impurities in a quantum many-body system. 



Recent developments in experimental quantum gas re- 
search focus on single particle control in a many-body 
system for detection and engineering of strongly corre- 
lated quantum states, including high-resolution imaging 
of strongly correlated bosonic systems in optical lattices 
(Tj-IH or the preparation of a mesoscopic number of degen- 
erate Fermions While impurities have been studied 
in balanced or imbalanced mixtures [5], extending this 
approach to single or few, individually controllable impu- 
rities in a quantum gas grants access to a huge number of 
novel applications proposed. In the direction of quantum 
information processing, atomtronics applications are en- 
visioned with single atoms acting as switches for a macro- 
scopic system in an atomtronics circuit [6] ; two impurity 
atoms immersed in a quantum gas can entangle by an ef- 
fective long-range interaction mediated by the gas 0, 
or individual qubits can be cooled preserving internal 
state coherence |8|, l9| . In the field of condensed matter 
simulation, in Bose gases strongly coupled Frohlich-type 
polarons have been predicted to form [l^. Our system 
paves the way for experimental studies of multi-polaron 
systems of one to ten impurities not only in the weak 
and intermediate coupling regime [ll^ , but also for strong 
coupling [lij. Adding impurities, and hence polarons, 
one-by-one would allow to experimentally track the tran- 
sition even to the many-body regime, and moreover yield 
information about spatial cluster formation . Applied 
to fermi systems, well localized single spins might allow 



for realizing model systems to study Kondo physics [14 1. 
Finally, fundamental questions of quantum physics can 
be addressed with unprecedented precision, as single im- 
purities can act as local and non-destructive probes of 
strongly correlated quantum many-body states [15]. Sin- 
gle atom sensitivity allows to observe single events of 
molecule formation yielding valuable insight into interac- 
tion properties, as we show in this paper, similar to work 
using single trapped ions [l6|. Further, adding single 
impurities one-by-one to an inititally integrable system, 
such as a quasi ID Bose gas [l7|, allows to controllably in- 
duce thermalization of a non-equilibrium quantum state. 

In addition to tight control over individual impurities, 
it is a well controlled interaction at ultracold tempera- 
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FIG. 1. Schematic setup of the experiment, (a) Rb atoms 
are precooled in a magnetic trap and then transferred into 
a crossed dipole trap to a position in close vicinity of a 
single-atom high gradient magneto-optical trap for Cs. Sub- 
sequently, both species are loaded into separated sites of a 
IDl lattice, (b) By adiabatically ramping down the lattice, 
single Cs atoms are immersed in the ultracold Rb gas stored 
in the crossed dipole trap. 



tures between quantum gas and impurity which lies at 
the heart of all these applications in order to preserve 
the coherence of the sub-systems. In recent realizations 
of single trapped ions immersed in a BEG [l8|, mK ion 
temperatures currently obstruct a coherent time evolu- 
tion of the hybrid system. 

Here, we present the deterministic immersion of single 
or few neutral Gs atoms in an ultracold Rb gas. We ex- 
plicitly demonstrate the time-resolved sympathetic cool- 
ing of a well-defined number of impurity atoms down to 
single impurities to the temperature of the many-body 
system, governed by elastic 5-wave colHsions as required 
in all scenarios above. As a first application exploiting 
the properties our hybrid system, we investigate three- 
body collisions event-by-event with single atom resolu- 
tion, yielding a precise value of the three-body decay co- 
efficient in Rb-Rb-Gs molecule formation. 

We use an ultracold gas above the condensation thresh- 
old at a high phase-space density of about 0.2. This sup- 
ports a clear interpretation of the observed interaction 
effects in the absence of quantum effects. However, the 
quantum regime will be at the focus of future experi- 
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merits, and it can be easily accessed by cooling the Rb 
gas a bit further below the critical temperature. 

The insertion of single ^^^Cs atoms relies on methods 
described in Ref. [l9|. An ultracold ^^Rb gas is produced 
and stored in a magnetic-field insensitive state (F = 1, 
mp = 0) in a crossed dipole trap, where F and mp are 
the total angular momentum and its projection onto the 
quantization axis, respectively. Single Cs atoms are cap- 
tured in close vicinity in a high gradient magneto-optical 
trap (MOT). Subsequently, both species are loaded into 
a one-dimensional (ID) optical lattice at separated lat- 
tice sites (Fig. [T]), and the Cs-MOT is switched off. Here, 
the Cs atom has a temperature of ~ 30 pK. Finally, the 
lattice is adiabatically removed while keeping an over- 
all trapping potential by the crossed dipole trap, forcing 
both species to interact in this common potential. Adi- 
abatic expansion upon removing the lattice cools the Cs 
atom to a temperature of ^ 5 pK in the crossed trap, and 
the temperature of the Rb gas is adjusted to either 250 
nK or 700 nK. The ultracold gas is detected by absorp- 
tion imaging is employed; impurity atoms are counted by 
monitoring their fluorescence after recapture in the MOT 
using a single-photon counter. The recapture probability 
for an atom present in the optical trap is close to unity, 
so that we assume that the fraction of atoms recaptured 
is identical to the survival probability. 

In balanced Rb-Cs mixtures, a large repulsive inter- 
species interaction was observed, which together with 
measured Feshbach resonances [20| led to a precise knowl- 
edge of the two-body molecular interaction potential (21I - 
|23|. A rapid thermalization of the immersed impurity is 
therefore expected. We study the thermalization of on 
average two "hot" Cs impurity atoms in F = 3, inter- 
acting via s-wave collisions with a large "cold" Rb cloud 
with A^Rb ^ 12000 in F = 1, mp = 0. As the lattice is 
removed, the impurity gets in contact with the Rb gas, 
and the sample is stored in the crossed dipole trap with 
a potential depth of about 30 pK (70 pK) for Rb (Cs). 

After a variable interaction time ti, Rb is pushed out 
of the trap with a resonant light pulse. We have verified 
that this light pulse does not influence Cs. For each in- 
terspecies interaction time ti, the temperature of Cs is 
measured with the release-recapture method 0, [lH . 

For this, the trap is switched off diabatically for a re- 
lease time At, in which Cs is allowed to expand freely. 
The survival probability after At depends on the kinetic 
energy and thus on the temperature: Hot atoms leave 
the recapture volume faster than cold atoms. The limit 
of survival is given by gravity which leads to a loss of 
atoms even at T = in a time of At\ ^2.5 ms for our 
trap parameters, see inset of Fig.[2fa). For atoms at finite 
temperatures exhibiting an energy distribution, this de- 
cay is smeared out towards smaller survival probabilities 
at shorter At. 

The inset of Fig. [2fa) shows two examples of release- 
recapture measurements for two different interspecies in- 
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FIG. 2. (a) Cooling of single Cs atoms interacting with an 
ultracold Rb gas. Within about 50 ms, Cs is cooled to the 
temperature of the Rb gas, the solid line shows the cooling 
according to our model. The horizontal line indicates the 
Rb temperature TRb ~ 250 nK, vertical dashed line indicates 
the time when the optical lattice is fully extinguished. Inset: 
Release-recapture measurement for Cs without (open circles) 
and with 117 ms interaction time with Rb (closed circles). 
Solid lines indicate the corresponding simulation. The dashed 
line shows the simplified T==0 recapture curve, see text, (b) 
Lifetime of Cs during cooling, (c) Normalized density distri- 
bution of Rb (light shaded) and Cs before sympathetic cooling 
(dashed) and after sympathetic cooling (shaded). 

teraction times tj. The survival probability starts at 
about 80 %, limited by losses during the loading pro- 
cedure pJi]. With increasing At the survival probability 
decays in a characteristic way. Cs atoms not interact- 
ing with Rb show a release half-life time of ~ 1 ms, 
whereas Cs atoms interacting for 117 ms with Rb feature 
an almost doubled release half-life time of ^ 2 ms, which 
indicates cooling. 

To quantitatively estimate the temperature of the im- 
purity atoms, the release-recapture experiment is mod- 
eled by a numerical simulation. The inset of Fig. [2fa) 
shows the corresponding simulation that best fits the ex- 
perimental data. For every interaction time the main 
graph of Fig. [2fa) shows the extracted temperature. 
From the initial temperature of T ^ 4.8 pK without in- 
teraction, Cs is cooled within about 50 ms to the tem- 
perature of the Rb gas of T ^ 250 nK, as determined 
independently by time-of-flight velocimetry. For even 
longer interaction times, the temperatures of Cs and Rb 
agree within the uncertainty, which we interpret as ther- 
mal equilibrium between the two sub-systems, forming 
an ideal starting point for the realization of the scenarios 
listed above. Our system moreover also yields insight into 
the non-equilibrium cooling dynamics, where in principle 
we have access to the full energy distribution at all times. 

In order to model the sympathetic cooling over 
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the full range, the interspecies density overlap n = 
^N^'^N^'' I ^Rb^Cs<^^^ and the mean relative velocity v = 



sku ( jRb_ _^ Tc^ \ ^j.^ assumed to be time-dependent. 

TT ymRb mcs J ^ 

The interspecies scattering cross section dRbCs can be 
calculated from the thermalization by 
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AT = CTRbCs^^ 



(1) 



where ^ = 0.96 is the reduction factor due to the 
mass difference [l^. Fig. [2ja) shows the result of 
solving Eq. [T] numerically for an energy independent 
scattering cross section crRbCs, that best fits the exper- 
iment. The shape clearly deviates from a purely expo- 
nential cooling as usually observed in the thermalization 
of bulk gases [ij, Hsl, which is a consequence of the 
negligible perturbation of the Rb gas by the impurity. 
The corresponding effective 5-wave scattering length is 

|0^RbCs| = \/crRbCs/47r ^ 450 ttQ. 

The s-wave scattering lengths for Rb-Cs collisions 
are known to be as = (997±ll)ao and ax = 
(513.3 ± 2.2) gp for the singlet and triplet potentials, re- 
spectively [23| . Since we only control the hyperfine state 
of Cs and Rb is in the F = 1, mp = 0-state, |aRbCs| 
is comprised by contributions of scattering channels of 
the singlet and triplet potentials, so that one would ex- 
pect ax < I^RbCsl < cis- The discrepancy of our result 
can be explained by two issues particular to our exper- 
iment. Firstly, we load Cs off-center into the trap so 
that initially the impurity atom could follow a trajectory 
which has a reduced overlap with the Rb gas. Secondly, 
the adiabatic lowering of the lattice causes Cs atoms to 
leave their initial site at non-zero lattice depth, similar 
to standard evaporation. We indeed observe cooling of 
Cs already a few ms before completely extinguishing the 
lattice, indicated by the vertical dashed line in Fig. Efa). 
These two effects lead to a decreased effective interaction 
time with respect to ti. Therefore, the derived value for 
|<^RbCs| should be regarded as a lower limit and the sym- 
pathetic cooling can be understood in the framework of 
the well-known two-body collision properties 2ll-[23| . 

The dynamics of our thermalization experiment are 
remarkable: During sympathetic cooling, the survival 
probability of impurity atoms decays by a few percent 
only, as illustrated in Fig. [2fb). Since we precisely know 
the initial number of impurity atoms, we observe the in- 
teraction of individual particles, rather than a statistical 
ensemble of a bulk gas in imbalanced mixtures [HI, [27| . 
This could be employed, for instance, for the preparation 
of a mesoscopic number of ultracold or degenerate par- 
ticles. The heating of the Rb gas, caused by the sympa- 
thetic cooling of a single impurity Cs atom, is estimated 
to be below one nK and can be neglected. Furthermore, 
the atom number of the Rb gas is approximately con- 
stant (discussed below), thus the properties of the Rb 



gas remain unaffected and the many-body system is not 
perturbed by the impurity. Cs, in contrast, is cooled 
by more than one order of magnitude in temperature, 
causing the trapping volume to shrink significantly. In 
Fig. EJc), the normalized density distributions of both 
species before and after thermalization are plotted, il- 
lustrating the immersion of the impurity atom into the 
Rb gas. Thus, our results represent the realization of a 
non-perturbative temperature (and density) probe for a 
many-body system, which could be extended to a local 
and coherent probe. 

While interspecies elastic two-body collisions, that de- 
fine t he prop erties of our hybrid system, are well under- 
stood 
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-23|, interspecies inelastic three-body collisions, 
limiting the lifetime, are more challenging to handle theo- 
retically and experimentally. The combination of Rb and 
Cs is particularly involved, as this mixture shows near- 
resonant scattering properties. Both the Cs-Cs and the 
Rb-Cs interaction potential feature weakly bound dimer 
states, which cause a very high three-body recombination 
rate in both loss channels. Three-body recombination 
then causes significant heating, leading to non- negligible 
evaporation loss from the trap (" ant i- evaporation" |28|), 
which in general affects each species to a different level, 
due to species-dependent trap depths. These effects lead 
to complicated dynamics, i.e. time-dependent changes in 
density and temperature of each species, impeding a clear 
analysis of the three-bod y p rocesses in (nearly) balanced 

Additionally, the analysis 
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of losses in bulk gases is further hindered by avalanche 
effects jaU. 

Our hybrid system represents a novel approach that is 
not affected by these effects. We can employ single impu- 
rities as individual probes that allow observing inelastic 
loss mechanisms atom- by- atom and event-by-event. 

Preparing both species in their respective absolute 
groundstates (Rb: F = l,mF = 1, Cs: F = 3,mF = 3) 
allows for inter- and intraspecies three-body losses, only. 
After preparation, both species are stored in a crossed 
dipole trap at a temperature of T ^ 700 nK, with a 
potential depth of 54 pK (96 pK) for Rb (Cs). Due 
to the wavelength of the dipole trap laser, the den- 
sity distributions are approximately the same for both 
species, and the differential gravitational sag is negligi- 
ble for the following experiments. The Rb peak density 
is 1.2x10^^ cm~^ and the three-body loss coefficient is 
Ks = 3 X 10~^^ cm^s~^, so that inelastic intraspecies Rb 
three-body collisions can be neglected on the relevant 
time scale, see Fig. [3jb). Furthermore, the data shows 
that the lifetime of the Rb gas is not affected by the im- 
purity atoms, as only few atoms can be lost causing a 
relative Rb loss of ^ 10~^. 

For the impurity atoms. Fig. [3fa) shows the lifetime 
averaged for 1 - 4 Cs atoms immersed in the Rb gas, well 
described by a simple exponential decay with time con- 
stant Tt = (732 ± 47) ms. The atomic resolution of the 
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impurity atom number allows us to post-select the data 
for a precisely defined atom number. The corresponding 
data extracted from the averaged measurement is shown 
in Fig. [3fc). For a single impurity atom, the only remain- 
ing loss channel is a Rb-Rb-Cs collision. Accordingly, for 
two impurity atoms also Cs-Cs-Rb collisions are possible, 
whereas for three and four impurity atoms all loss chan- 
nels are open. With the same decay as obtained from the 
average data, the lifetimes of a single, two, three and four 
atoms are well described, suggesting that Cs-Cs-Rb and 
Cs-Cs-Cs three-body collisions can be ruled out here. 

This can be further verified by analyzing the statistics 
of the losses, relevant for future studies of two-body in- 
teractions. For the case of exactly two impurity atoms 
immersed in the Rb gas, the respective probabilities pi, of 
losing i = 0, 1, 2 Cs atoms are determined for each inter- 
action time ti, corresponding to a certain total survival 
probability p. Assuming independent impurities and thus 
the same, independent total survival probability p for 
each Cs atom, these probabilities are given by po = p^, 
Pi = 2p{l — p) and p2 = (1 — — p)- In Fig- IH this 
expectation is plotted together with the corresponding 
relative occurrences obtained from the experiment, show- 
ing good agreement. The impurity atoms are therefore 
independently interacting with the ultracold Rb gas and 
hence the observed decay is due to Rb-Rb-Cs three-body 
recombination. We stress that the ability of studying ev- 
ery single interaction event with single-atom resolution 
is the key for the unambiguous interpretation of inelastic 
collisions. 

Overall, the decay of survival probability of Cs atoms 
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FIG. 3. (a) Decay of on average 1-4 impurity atoms. The 
solid line is an exponential decay, (b) Rb atom number for the 
same time interval. The solid line is the theoretical lifetime 
due to intraspecies three-body decay and collisions with the 
background gas. (c) Lifetime of a precisely defined number of 
impurity atoms, post-selected from (a). Solid lines show the 
same fit as in (a). The insets illustrate the respective open 
loss channels. 
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FIG. 4. Three-body losses for exactly two impurity atoms, 
observed event- by- event, (a) Fluorescence traces showing all 
three possible loss events and the according allowed loss chan- 
nels, (b) Corresponding relative occurrences for all possible 
case: two atom loss (a), one atom loss (t) and no atom loss 
(•). The solid lines show the expectation, for details see text. 



due to three-body collisions with Rb can be described by 
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Ncs 



(2) 



where (n^^) is the mean squared Rb density and Nqs 
the Cs atom number. As described above, the mean 
squared Rb density {n'^^) is essentially constant and 
Eq. (|2j) is simply solved by N{t) = Nq ex.p {-t / n) = 
A^oexp (— L3(n^^)t). We obtain the three-body loss co- 
efficient Ls = (5 ± 2) X IQ-^^cm^s-^ [111, in very good 
agreement with a recent theoretical calculation [33| . Our 
approach allows therefore a precise determination of L3, 



so far impossible with balanced mixtures [22 



In conclusion, our experiment demonstrates the cre- 
ation of a novel hybrid system, comprised of an atomic 
impurity interacting with a many-body system via 5- wave 
collisions. Recently found Feshbach resonances [2o| will 
allow us tuning the interaction of the two subsystems in 
future experiments, while the use of species selective po- 
tentials [3J] will improve localization of impurities with- 
out perturbing the gas. Together with site-resolved de- 
tection beyond the diffraction limit [35|, this provides a 
good starting point for the realization of the scenarios 
mentioned above. 
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